nmds <-metaMDS(abundance_matrix, distance ="bray", k =2)
Wisconsin double standardization
Run 0 stress 0.08532757
Run 1 stress 0.1200854
Run 2 stress 0.08703641
Run 3 stress 0.09333725
Run 4 stress 0.08902399
Run 5 stress 0.1081445
Run 6 stress 0.08492605
... New best solution
... Procrustes: rmse 0.01883282 max resid 0.2615482
Run 7 stress 0.08791272
Run 8 stress 0.1139224
Run 9 stress 0.0854949
Run 10 stress 0.1158335
Run 11 stress 0.119714
Run 12 stress 0.09171769
Run 13 stress 0.08640987
Run 14 stress 0.08262548
... New best solution
... Procrustes: rmse 0.008558615 max resid 0.06322369
Run 15 stress 0.08490224
Run 16 stress 0.08386775
Run 17 stress 0.1077224
Run 18 stress 0.08189793
... New best solution
... Procrustes: rmse 0.008302257 max resid 0.0633948
Run 19 stress 0.0937349
Run 20 stress 0.08654519
*** Best solution was not repeated -- monoMDS stopping criteria:
10: no. of iterations >= maxit
8: stress ratio > sratmax
2: scale factor of the gradient < sfgrmin